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Abstract 

For the simplest minisuperspace model based on a homogeneous, isotropic metric 
and a minimally coupled scalar field we derive analytic expressions for the caustic 
which separates Euklidean and Minkowskian region and its breakdown value 0* . This 
value represents the prediction of the no-boundary wave function for the scalar field 
at the beginning of inflation. We use our results to search for inflationary models 
which can render the no-boundary wave function consistent with the requirement of 
a sufficiently long inflationary period. 



Quantum cosmology U can shade new light on the problem of initial conditions in the 
universe. To be predictive boundary conditions have to be imposed on the wave function 
of the universe. What quantum cosmology can do is therefore to replace the procedure of 
imposing specific boundary values by a "theory of boundary conditions" . Several proposals 
for such a theory have been made. The most prominent ones are the no-boundary proposal 
due to Hartle and Hawking and the tunneling proposal due to Linde and Vilenkin [|J. 
Clearly, all these proposals are speculative and have to be confronted with reality i. e. with 
the main characteristics of our universe. 

In the simplest minisuperspace model which is based on a homogeneous, isotropic met- 
ric with scale factor a and a single, minimally coupled scalar field — 0W both the 
no-boundary and the tunneling wave function predict a correlation between a, and their 
conjugate momenta if a is sufficiently large. Assuming decoherence this correlation corre- 
sponds to a classical universe in a de-Sitter phase. They are therefore in agreement with 
the classical behaviour of the universe we live in and the inflationary paradigm of cosmol- 
ogy. In addition, a probability distribution P(0) for the field at the beginning of the 
inflationary period can be extracted from the wave function. It can be used to compute a 
conditional probability P su g for sufficient inflation : 



Here su g is the minimum value of which guarantees the A^ min ~ 60 e-folds of inflation 
needed to solve the problems of standard cosmology. Since the wave function has been 
determined in a semiclassical approximation one should cut off the integral at the value 02 
corresponding to Planck energy. <f>i is the minimal which leads to an expanding universe. 
For the tunneling wave function the probability distribution 



prefers large values of and therefore generically predicts sufficient inflation. However, the 
opposite is true for the no-boundary result 



Since conventionally the potential V contains a tiny coupling in order to keep the density 
fluctuations small enough Pnb is so strongly peaked on low values of that it essentially 
predicts the lower cutoff value 0i which is much smaller than su ff. Clearly, this argu- 
mentation can be doubted since the behaviour of P(0) beyond Planck energy where the 
above expressions are not valid may change the picture completely. It has, however, been 
argued that the one loop contribution to the distribution P((f>) can suppress the constant 
tail at large PJ. This can render P(4>) normalizable and serve as a justification for the 
semiclassical approximation. Here we assume that such a suppression indeed takes place 
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which, in addition, leaves the low energy behaviour of the probability distribution - as 
described by the above equations - essentially unchanged. Furthermore, we assume that 
we stay in a regime where the back reaction of quantum fluctuations can be neglected and 
the field follows a classical slow roll trajectory. Then P su g should indeed be large for a 
successful inflationary model and its computation provides a check for a certain boundary 
proposal. 

Grishchuk and Rozhansky || have pointed out that the existence of a caustic in the 
a-0-plane is a necessary condition for a transition from the Euklidean region at small a 
to the classical Minkowskian region at large a. For small this caustic breaks down so 
that below a certain value 0* > 0i no classical universe develops. Therefore, in eq. ([!]) 0i 
should be replaced by 0*. Then the no-boundary wave function predicts 0* as the value of 
at the beginning of inflation. This value strongly depends on the potential and it has to 
be clarified whether 0* can be larger than su g for some choices. 

In ref. [[J 0* has been determined numerically for the potentials V = A0 2 and V = A0 4 . 
In both cases it was found to be smaller than su ff . However, it was left as an open problem 
whether or not a potential with 0* > su g- exists for the no-boundary wave function. This 
is the question we are going to address in this paper. As a first step we will derive an 
equation for the caustic and analytic expressions for 0*. These expressions allow us to 
discuss the general dependence of 0* on the potential V(4>) and to analyze several classes 
of potentials in detail. 

The Euklidean equations of motion for a and which have to be solved in order to 
determine the semiclassical wave function are 
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<t> + 3-0 = --r- 
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According to the no-boundary proposal they have to be integrated by using the initial 
conditions 

a(r = 0) = , 0( r = 0) = , 0(r = 0) = O . 

To illustrate the qualitative behaviour in fig. [I] we have plotted the trajectories for a number 
of starting values O using the potential [] V = A0 4 . The caustic consists of the points of 
maximal a and in this specific example it breaks down at 0* ~ 3. For large values of it 
behaves asymptotically like 

a 2 V = 1 (5) 

(dashed line in fig. P as it can be expected from the second eq. (§). Near the breakpoint, 
however, it deviates substantially from this asymptotic expression. As a first step we will 
now improve eq. ([|) such that it also describes the caustic for small 0. 



L The value of </>* does not depend on the coupling A which we have set to 1 for the numerical integration. 
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Figure 1: Trajectories for the potential V = 4 , approximate caustic a 2 V = 1 (dashed line) 
and improved caustic (dotted line). 

Obviously, informations about the time dependency of a and are not needed for this. 
Instead, we can look at the equation 



(1 - a 2 V)<\>" + 



-(3 - 4o 2 V) - a(2 - 3a 2 VU' 2 
a 



1-aV 2 dV 



-n (6) 



for = 0(a) (0' = dcf)/da) which can be derived from eq. Surprisingly, this equation 
has two exact solutions 

0' = ±- (7) 
a 

which are independent on the potential V. It turns out that the second of these solutions 
(0' = —1/a) asymptotically describes the "return trajectories" i. e. the trajectories after 
the turning point at maximal a. Since we are mainly interested in the behaviour near the 
approximate caustic a 2 V = 1 we linearize eq. (H) in e = a 2 V c — 1, V c = V(0 C ) to obtain 



de \ AJ \ 2y c J de 
We have introduced the variable y 
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which will be quite useful in the following. Eq. (§) can be integrated explicitly to give 



G(0 
G(0 




(10) 
(11) 



with 



e c (j/c 



(12) 

where we have imposed the boundary condition £ = 0ate = — 1. At e = e c (y c ) our solution 
shows a singularity £ — > oo which we identify with the turning point of the trajectory. From 
this observation we get for the caustic : 



1 + CciVc 



(13) 



In fig. |T] this curve is shown for our example V = 4 (dotted line). 

To determine 0* we have to distinguish two cases. Let us first assume that the potential 
in the range we are interested in is increasing rapidly. This is e. g. fulfilled for a polynomial 
potential like the one in fig. |l|. Then the first part of the trajectories is approximately 
constant in and because of the strong variation of the caustic with a the trajectories 
follow their asymptotic behaviour 0' = — 1/a even before they intersect the caustic for the 
second time. Then an expression for 0* can be computed by the following prescription. 
A return trajectory defined by 0' = — 1/a which starts at some point (a c , C ) on the 
caustic ( fL3|) intersects the approximate caustic a?V = 1 at (a,(f>). Using the eqs. (|5]), (|7|), 
( JT~5| ) and ( |T2"D we find the relation 
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1 + e c (y c 
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Now we minimize as a function of C to identify the lowest return trajectory (cf. fig. 1). 
The value of 0(0 C ) at the minimum can be identified with 0*. The explicit computation 
results in conditions on y c = y(0 c ) and 0* : 



1 

Vc 



1 de c / dy 
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(y c ) y'(y c ) = i 
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1 + e c (y c ) 



(15) 
(16) 



From eq. (|I2| ) it can be seen that e c (y = 1) = e 1 and de c /dy(y = 1) = — e 1 . Therefore, if 

y\y) = 2(l -\y^{y)\ (17) 
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is ~ 0(1) at y = 1 we can consistently neglect the second term in eq. (0). We will see that 
this is indeed justified for a number of potentials. The above equations then simplify to 

2V 

Vc = -^(0c) = l (18) 
= k^=. (19) 

with 

k = Vl + e- 1 ~ 1.17 . 

These equations now allow an explicit calculation of 0* for a given sufficiently increasing 
potential. Since the constant k in eq. fll9| ) is close to unity C and 0* are of the same order 
for almost all potentials 0. An order of magnitude estimate for 0* is therefore given by the 
solution to the simple equation 

2y , . 

V = y? = 1 ■ ( 2 o) 



For potentials with a very flat region e. g. caused by a maximum or an asymptotic 
behaviour V — > const the situation is different. Fig. |^ shows the example V — (1 — 
exp(— 4>/f)) 2 for / = 0.1. The trajectories increase substantially from a = up to the 
caustic which is nearly vertical. The return up to the second intersection with the caustic 
is too short to be described by the asymptotic equation 0' = —1/a. Instead, we follow a 



trajectory starting out at yo 
using the equations 



ywo) 



and intersecting a 2 V — 1 at y c — y((p c ) by 
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They can be derived from eq. (||) and from eq. (|10D by expanding G(£) up to the first order 
in l/£. These approximate first integrals for small a and the region around the caustic 
can be "glued" together. Their integration gives the second intersection value with the 
1 as a function of 0o- As before the minimum of this function should be 
0*. Following this prescription we arrive at a set of equations 
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(21) 



l 



y'(yo) 



2 An. exception to this is given by an exponential potential which cancels the exponentials in the numer- 
ators of eq. (|l9|) to a high accuracy. We have discussed this special case in an example below. 
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Figure 2: Trajectories for the potential V = (1 — exp(— 4>/ f)) 2 for / = 0.1, approximate 
caustic a 2 V = 1 (dashed line) and improved caustic (dotted line). 

where we have used that e c (y) ~ 1/y 2 since y ^> 1 in a flat region of the potential. 

For the two types of potentials - the increasing "chaotic" type of potentials and the 
type with a flat region - the eqs. (|T5|), ([16]) and eq. fl2~T|) represent the general answer for 
the value of </>*. Clearly, we have not proven that our approximations are justified for all 
potentials. The results should therefore be used as a guide to find promising potentials 
whose properties then have to be verified numerically. In this way large classes of potentials 
can be analyzed without a numerical case by case study. In fact, for all potentials we have 
checked so far our analytic formulae gave the correct answer. 

We start by discussing the consequences of the previous results on the first type of 
potentials. The number of e-folds in the slow roll approximation is given by 

iV e _ folds = 3 (*' y{<t>) d<j) (22) 



where 0/ is the end point of the slow roll regime determined by the breakdown of the 
conditions 

1 V" 

*/»-, -^«18- (23) 

Suppose we have found a potential V of the first type which generates the desired number 
of e-folds. Since y does not increase "irregularly" it is plausible that the slow roll condi- 
tions (|23|) are still fulfilled at C . Consequently, our approximation leading to (0) and ( |19D 
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is justified and we have y c ~ 1 at C . This leads to the rough approximation 

iV e _ folds < 30* . (24) 

It shows that for such potentials one will need a quite large value 0* to generate the needed 
iVmin ~ 60 e- folds. From y(0) < 1 for ~ 0* we conclude that 

> p 2(0 2 - 01 ) 

for any 02, 0i ~ 0*. In particular we have 

V{4>*) * Ae 2 ** ~ Ae 2 ^ ~ Ae 40 (25) 



where A is the coupling in front of V. Therefore a potential which can provide a large 
enough 0* either exceeds Planck energy at this point - which makes the predicted 0* 
unreliable - or needs a very small coupling A. In the second case quantum fluctuations 
in de Sitter space are generically too small to generate a sufficiently large initial density 
spectrum. Then other sources of such density perturbations have to be considered. 

Examples are provided by the standard chaotic potentials V = A0 n , n — 2, 4, 6, • • •. It 
is easy to show that y'(y = 1) = 2/n so that the eqs. (|T^) and (|TJ) can be applied. We get 

0c = \ ■ (26) 



Writing 0* = C + 5 for sufficiently large n we can approximate y V((j)*) ~ JV((f) c )(e(l 
5/n)) s to obtain an explicit solution of eq. ([TJD : 



+ Vl - k- 1 y/n — — + 0.38v^ . (27) 



2 2 

By numerical integration of eq. (|^) we have checked up to n = 70 that this formula indeed 
results in reasonable values for 0*. The number of e- folds is given by 

iVc-foids ~-4>l~-n + 1.14v^ . 
n 4 

In accordance with the general discussion a very large value of 0* and correspondingly of 
n is needed to generate sufficient inflation. In fact this is possible for V = A0 n with n ~ 70 
but an enormously small A is needed to keep V(<f)*) below Planck energy. 
Let us next discuss the potential 

V- = A 4 (l-exp(^ (28) 

which (for / = 1) can be derived from higher order gravity theories 0. We concentrate on 
the interesting range / ~ 1 where the exponentials in eq. (|1^) are canceled. In that case 
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0* may be substantially larger than C and the above general argument is somewhat weak. 
From the eqs. flT|), ([Ig) and @ we obtain forO<d = /- l«l: 

0c ^ In 3 

a 

, 1 In k , . 

0, ~ ln- + — (29) 
a a 

iV(0,) ~ 30, . 

These formulae have been confirmed by numerical integration. For d ~ 0.01 a sufficient 
number of e-folds may be generated, however, in accordance with our general conclusion a 
small coupling A 4 ~ e -40 is needed to keep V(0*) below the Planck scale. 



ln + 4 - 20 / ■ ( 3 °) 




Now we turn to the second type of potentials with a flat region. An interesting example 
is provided by the above potential (|2"8"| ) in the negative direction where it approximates 
a constant value. Quantum cosmology of this potential has been studied in a number of 
papers including 0] . Only for a value of 0* in the asymptotically flat region we can hope for 
sufficient inflation. Such values are obtained for / <1 where we can solve the system ( pT|) 
explicitly. We find : 

1 

/ " ~ 1.45/ 2 

This has been confirmed numerically. Given the formula 

A^e-folds( 

we see that independent on / only a few e-folds are generated starting at 0*. 

Another way to obtain a flat region is to consider a potential with a maximum at some 
value = m . A sufficiently long period of inflation will be produced for these potentials 
if starts its evolution in a small range around the maximum. Basically, we have to ask 
two kinds of questions : Can the shape of the potential be such that 0* falls into this range 
and is there a fine tuning involved in arranging such a situation ? 

Before we address these questions in general let us be more specific and discuss a 
prominent example for these potentials : The natural inflation potential V = A 4 (l — 
cos(20//)) which can be interpreted as the effective potential of a pseudo-goldstone boson. 
Its inflationary properties have been studied extensively in ref . || . The quadratic maximum 
of this potential is still "curved" enough to apply the eqs. fll8D and (0). From 

(24>c\ f 2 -i ,,,, 

cos It] = —i (31) 

we see that C and correspondingly 0* approach the maximum for decreasing /. The latter 
is found from a numerical solution of eq. and is in agreement with the results from 
numerical integration of (01). At / ~ 0.6 the value 0* reaches the maximum. In principle, 
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one can therefore hope that a value of / close to 0.6 produces sufficient inflation. However, 
taking into account that 

, - _ „ exp (-^j (32) 

this seems almost impossible : For / ~ 0.6 the starting point of the field <fi has to be 
at the maximum to such an enormous precision that one has to fine tune / in order to 
obtain a 0* in this range. Even if such a tuning is accepted one would have to worry about 
the fluctuations of the higher, inhomogeneous modes and their back reaction on the field 
0. This argument does not depend on the specific choice of the potential but just on the 
quadratic behaviour around the maximum. Our conclusion holds for all potentials with 
such a quadratic maximum provided that the slow roll near the maximum is the main 
source of inflation. 

We can generalize this analysis to higher order maxima. If we assume that the relevant 
phenomena appear near the maximum we can expand the potential around m up to the 
first nontrivial order In : 

V = A 4 (1-A 2 ") + 0(A 2 ™ +1 ) , A = ^y^. (33) 

For the number of e-folds we get 

iV e -foids(A*) ~ 2n *f_ 1} A l_ 2 for n > 1 . (34) 

In contrast to the n — 1 case the number of e-folds does not show a logarithmic dependence 
on the initial value so that fine tuning problems are irrelevant. The equations ([H]) can be 
solved with the ansatz y = d n /f n ^ n ~ 1 \ y c = Cn/ f n /( n ~ l ) leading to 

A, ~ (nc n y^f^ (35) 

3 In — 2 

iV e _ folds (A,) ~ — -(nc n )— . (36) 

2n(n — 1) 

As for the exponential potential the number of e-folds turns out to be independent on /. 
The numerical constants c n are small, e. g. C2 — 2.6, C3 ~ 5.0, showing that A r e -f ids(A*) 
remains below the needed value. For n = 2, 3 this result has been verified by numerical 
integration. 



In conclusion, we have determined an analytic equation for the caustic separating the 
Euklidean and Minkowskian region in the simplest minisuperspace model. Moreover, we 
derived formulae for the breakdown value 0* of the caustic which is the value of the scalar 
field predicted by the no-boundary wave function. Our results give some insight in the 
structure of minisuperspace models and its dependency on the scalar field potential. They 
allow to analyze inflationary properties of these models without relying on a numerical 
case by case study. We have applied these results to a wide range of models with a single 
scalar field. Unfortunately, for potentials with an extraordinary flat region, e. g. caused 
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by a maximum or an asymptotic behaviour V — > const we could not find an example 
which leads to a sufficiently long period of inflation. For potentials with a maximum and 
an inflationary phase which mostly arises due to the slow roll near the maximum this 
behaviour seems to be general. 

On the other hand, some potentials of "chaotic" type like V = \<p n with n ~ 70 can 
cause sufficient inflation. However, to stay within the limits of the semiclassical approxima- 
tion the coupling A generically has to be very small. A positive answer within these limits 
therefore asks for an alternative explanation of primordial density perturbations. From 
our results, one might also argue that an understanding beyond the semiclassical limit is 
necessary to decide about the existence of a "successful" model. On the other hand, we 
have concentrated on the simplest class of models with a single scalar field which in the 
inflationary phase behaves essentially classical. Non negligible quantum fluctuations, in 
particular in the context of multidimensional models, might, however, change the picture. 
From the viewpoint of particle physics models with more than one scalar field are the more 
realistic anyway. Generalizations of the results presented in this paper might be helpful in 
analyzing the quantum cosmology of such models. 
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